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TORSION AND GROUND STATE MAXIMA: 
CLOSE BUT NOT THE SAME 

B. A. BENSON, R. S. LAUGESEN, M. MINION, B. A. SIUDEJA 


Abstract. Could the location of the maximum point for a positive solution of 
a semilinear Poisson equation on a convex domain be independent of the form 
of the nonlinearity? Cima and Derrick found certain evidence for this surprising 
conjecture. 

We construct counterexamples on the half-disk, by working with the torsion func¬ 
tion and first Dirichlet eigenfunction. On an isosceles right triangle the conjecture 
fails again. Yet the conjecture has merit, since the maxima of the torsion func¬ 
tion and eigenfunction are unexpectedly close together. It is an open problem to 
quantify this closeness in terms of the domain and the nonlinearity. 


1. Introduction 

Suppose the Poisson equation 

{ —Am = f{u) in O, 
u = 0 on 

has a positive solution on the bounded convex plane domain hi. Here the nonlinearity 
/ is assumed to be Lipschitz and restoring, which means f{z) >0 when z > 0. Cima 
and Derrick [2, 3] have conjectured that the location of the maximum point of u is 
independent of the form of the nonlinearlity /. 

This conjecture sounds impossible, since the graph of the solution must vary with 
the nonlinearity. Numerical computations by Cima and co-authors give surprising 
support for the conjecture, though, and Figure 1 provides further food for thought 
by considering a triangular domain and plotting the level curves and maximum point 
for the choices f{z) = 1 and f{z) = Xz. The corresponding linear Poisson equations 
describe the torsion function and the ground state of the Laplacian (see below). Our 
solutions were computed numerically by the finite element method on a mesh with 
approximately 10® triangles. The maximum points for the two solutions in Figure 1 
appear to coincide, even though the level curves differ markedly near the boundary. 

We disprove the conjecture on a half-disk in Section 2, and again on the right 
isosceles triangle in Section 3. Interestingly, the conjecture is remarkably close to 
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Figure 1. Level curves and the maximum point on a triangular do¬ 
main, for solutions of two different Poisson type equations: the torsion 
function (left) and the hrst eigenfunction (right). 


being true in these counterexamples, with the maximum points occurring in almost 
but not quite the same location. We cannot explain this unexpected closeness. 

A fascinating open problem is to bound the difference in location of the maximum 
points of two semilinear Poisson equations in terms of the difference between their 
nonlinearity functions and geometric information on the shape of the domain. Also, 
note that for both the half-disk and right isosceles triangle, our results show that the 
maximum point of the torsion function lies to the left of the maximum for the ground 
state (when oriented as in Figure 1), which perhaps hints at a general principle for a 
class of convex domains. 


Notation. The torsion or landscape function is the unique solution of the Poisson 
equation 

—Am = 1 in ff, 

M = 0 on dQ. 

Here we have chosen f{z) = 1. Clearly u is positive inside the 
maximum principle. 

The Dirichlet ground state or first Dirichlet eigenfunction of the 
unique positive solution of 

—Am = Am in 
M = 0 on dQ, 

where A > 0 is the hrst eigenvalue of the Laplacian on the domain under Dirichlet 
boundary conditions. Here we have chosen f{z) = Xz. 


domain, by the 
Laplacian is the 
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Figure 2. The radial part of the ground state on the right half-disk: 
n(r, 0) = 


2. The half-disk 


The maximum points for the torsion function and ground state can lie so close 
together that one cannot distinguish them by the naked eye, as the following Propo¬ 
sition reveals. Yet the two points are not the same. 

Proposition 1. Take = {{x,y) ■. x > < 1} to he the right half-disk. 

On this domain the torsion function u attains its maximum at (0.48022,0) while the 
ground state v attains its maximum at (0.48051,0). Here the x-coordinates have been 
rounded to 5 decimal places. 

Proof, (i) The ground state is given in polar coordinates by 

v{r,9) = Ji{ji,ir) cos 9 

where Ji is the hrst Bessel function and ~ 3.831706 is its hrst positive zero. 
Clearly the maximum is attained on the x-axis, where 0 = 0, and the function is 
plotted along this line in Figure 2. By setting J((ji,ir) = 0 and solving, we hnd 
r = j[ — 0.48051, rounded to hve decimal places, where j[ ^ ~ 1.841184 is the 
hrst zero of J[. 

(ii) The torsion function is more complicated [5, Section 4.6.2], and is given by 



One verihes the Dirichlet boundary condition on the right half-disk by examining four 
cases: (i) n = 0 if x = 0 and 0 < |i/| < 1, (ii) u —)■ 0 as (x, y) —)■ (0, 0), (iii) n —)■ 0 as 
(x, y) —)■ (0, ±1), and (iv) u —)■ 0 as (x, y) —)■ (xi, yi) with xi > 0 and x^ yf = 1. 

To check u satishes the Poisson equation —An = 1, a lengthy direct calculation 
suffices. 

We claim u attains its maximum at a point on the horizontal axis. For this, hrst 
notice u is even about the x-axis by dehnition, meaning u{x,y) = u{x, —y). Hence 
the harmonic function Uy equals zero on the x-axis for 0 < x < 1. Further, < 0 
at points on the unit circle lying in the open hrst quadrant, since n > 0 in the right 
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half-disk and m = 0 on the boundary. Also, one can compute that Uy{x, y) approaches 
0 as (x, y) —)■ (0, 0) or (x, y) —)■ (1, 0) or (x, y) —)■ (0,1) from within the hrst quadrant 
of the unit disk. Lastly Uy vanishes on the y-axis for 0 < ?/ < 1 (since m = 0 there). 
Hence we conclude from the maximum principle that < 0 in the hrst quadrant of 
the unit disk, and so u attains its maximum somewhere on the x-axis. 

On the x-axis we have 


'(.( 1 , 0 ) = - 


27rx^ — 2x ^ + 2x + {2 + X ^ + x^) arctan 


2x 


x^ 


for 0 < a; < 1. Clearly m(0, 0) = u{l, 0) = 0, and 

rXx, 0) = [x + x^ - TTX^ + - 1) arctan • 




TTX-^ 


One can show by taking another derivative and applying elementary estimates that 
u{x, 0) is concave. Calculations show U 3 ;{x, 0) is positive at a; = 0.480219 and negative 
at a; = 0.480220, and so the maximum of u lies between these two points, that is, at 
X = 0.48022 to 5 decimal places. □ 


3. The right isosceles triangle 


Proposition 2. Take O = {{x,y) : 0 < a; < 1, ||/| < 1 —a;}, which is an isosceles right 
triangle. On this domain the torsion function u attains its maximum at (0.39168, 0) 
while the ground state v attains its maximum at (0.39183,0). Here the x-coordinates 
have been rounded to 5 decimal places. 


Proof, (i) Rotate the triangle by 45 degrees clockwise about the origin and scale up 
by a factor of tt/v^, then translate by 7r/2 to the right and upwards, so that the 
triangle becomes 

T = {{x,y) -. 0 < y < X < tt}. 

This new triangle has ground state 

v{x, y) = sin x sin 2y — sin 2x sin y = 2 sin x sin |/(cos y — cos a;) > 0 

with eigenvalue 1^ -|- 2^ = 5. One checks easily that n = 0 on the boundary of T, 
where i/ = 0 or a; = tt or y = a;. To hnd the maximum point, set = 0 and Vy = 
and deduce cos 2a; = cos a; cosy = cos2y. Therefore the maximum lies on the line of 
symmetry y = tt — a; of the triangle T. A little calculus shows that ^(a;, n — x) attains 
its maximum when x = arcsin(l/v^) + 7r/2. Hence the ground state of the original 
triangle attains its maximum at ((2/7r) arcsin(l/v^), O) = (0.39183,0) to 5 decimal 
places. 


(ii) The torsion function on the triangle T is 
u{x,y) 


= + 2 ^- 
n=l 


sinh nx sin ny — sin nx sinh ny 


2TTn^ sinh mr 

+ sinn(7r — a;) sinhn(7r — y) — sinhn(7r — a;) sinn(7r — y) 
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as we now explain. Observe that —Am = 1 because the inhnite series is a harmonic 
function, and m = 0 on the boundary of T by simple calculations with Fourier series 
when 0 < a; < TT, 1 / = 0, and when x = vr, 0 < i/ < tt; also m = 0 on the hypotenuse 
where y = x. 

The torsion function is known to attain its maximum somewhere on the line of 
symmetry y = ir — x, either by general symmetry results [2, 3] or else by arguing as 
in the proof of Proposition 1 part (ii). On that line of symmetry we evaluate 


u{x, n — x) 

“ — 2(1 — (—1)"') 

= —(x — ^- - [(—sinhnx — sinhnfvr — x)l 

^ vm^sinhuTT ^ ^ ^ u 

n=l 


sin nx. 


The series converges exponentially on each closed subinterval of (0, vr), and so we may 
differentiate term-by-term to hnd 


^U(X, TT — x) 
dx 


-2(x - 71/2) + 


00 2 

n 


TT^ - 2(1 - (-1)”) 


n=l 


Tim? sinhnvr 


( 1 ) 

coshnx -|- coshn(7r — x)] sinnx 


-f- [(—sinhnx — sinh?7,(7r — x)] cosnx k 


where once again the series converges exponentially on closed subintervals of (0,7r). 
The absolute value of the n-th term in series (1) is bounded by 


7r(e”^ + 

sinh(7T,7r) 


< 37r(e-"("-") + e-"") 


as we see by bounding the sin and cos terms with 1, adding the sinh and cosh terms 
having the same arguments, and using that sinh(n7r) > e"'^/3 for n > 1. Hence the 
inhnite series (1) is bounded term-by-term by dvr times the sum of two geometric 
series having ratios and e“*. 

The derivative of u along the line of symmetry is positive at x = 2.1860525 and 
negative at x = 2.1860530, as one hnds by evaluating the hrst 20 terms of the series 
in (1) and then estimating the remainder with the geometric series as above. Hence 
u has a local maximum at x = 2.186053 to 6 decimal places. This local maximum is 
a global maximum because y/u is concave (see [1, Example 1.1] or [4]). Translating 
to the left and downwards by 7r/2 and then scaling down by a factor of and 

rotating counterclockwise by 45 degrees, we hnd the torsion function on the original 
triangle has a maximum at 

X = -(2.186053 - 7r/2) = 0.39168 

TT 


to 5 decimal places. 


□ 
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4. Concluding remarks 

The counterexamples in this paper concern Poisson’s equation for f{z) = 1 and 
f{z) = Xz. One can hnd a whole family of counterexamples using f{z) = a+bz, where 
a > 0 and 0 < 6 < A. Note the maximum point depends on b but not a, as one checks 
by rescaling the solution u to u/a. To study this maximum point as b varies, one 
starts with the eigenfunctions of —A — b on the half-disk or right isosceles triangle and 
notes that the eigenfunctions are the same as for —A, just with eigenvalues shifted by 
b. The corresponding torsion function can be computed in terms of an eigenfunction 
expansion, and then the position of the maximum point can be carefully numerically 
located. We leave such investigations to the interested reader. 

Finally, while our counterexamples involve linear Poisson equations, our choices of 
/ could presumably be perturbed to obtain genuinely nonlinear counterexamples. 
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